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SUMMARY 


The  percentage  of  the  mss  flux  of  fine  particles  carried 
in  a  two  phase  fluid  that  impinge  on  the  vail  of  the  convergent 
section  of  a  nozzle  is  investigated  with  many  different  methods, 
to  reveal  the  qualitative  trends  of  different  physical  factors. 

If  the  particles  are  uniformly  distrl'outed  over  the  nozzle 
entrance  and  if  the  diameter  of  the  particles  is  less  than 
5  miciv^is,  the  particles  in  the  turbulent  core  outside  the 
displacement  area  of  the  turbulent  boundary  layer  on  the  wall  of 
of  the  rocket  chamber  will  pass  through.  The  percentage  of 
peirtlcle  mass  flux  that  Impinges  on  the  nozzle  wall  is  vexy  small. 

Qualitatively,  it  is  concluded  that: 

(a)  Most  of  the  impingement  takes  place  in  the  entrance  region, 
where  the  heat  transfer  condition  is  not  so  severe  as 

is  in  the  throat  region.  The  percentage  mass  flux  of 
particle  impingement  is  insensitive  to  nozzle  shape  near 
the  throat  but  depends  heavily  on  the  nozzle  geometry 
near  the  entrance. 

(b)  The  effect  of  moderate  changes  in  7  of  the  gas 
Is  negligible. 

(c)  Increase  in  particle  dispersed  density  causes  alight 
decrease  in  percentage  of  impingement,  but  the  absolute 
arx)unt  of  impingement  is  increased. 


For  the  sane  nozzle  geonetxy  and  particle  concentration, 
the  aoount  of  Implngesient  ia  directly  proportional  to 
the  particle  dianeter  squared. 

Increase  in  the  convergence  a’lgle  of  the  nozzle  causes 
em  increase  in  the  percentage  of  the  iiiQ>ingement. 
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1.  IirTRODUCTION 


In  many  modern  solid  propellant  rockets >  metallic  powders 
like  aluminxim  are  added  to  the  propellant.  The  metallic  oxides 
are  found  in  the  combustion  chamber  as  condensed  phases,  either 
as  liquid  droplets  or  as  solid  particles  of  micron  sizes.  The 
two  phase  fluid  expands  through  the  rocket  nozzle  with  con^llcated 
processes  of  beat  and  momentum  transfer  between  the  particles  and 
the  gar.  The  coi^densed  oxide  particles  (or  droplets)  are  much 
hotter  than  the  expanding  gas.  An  impingement  of  a  hot 
particle  on  the  nozzle  wall  will  not  only  decrease  the 
thrust -producing  outgoing  momentum,  but  will  also  increase  the 
heat  transfer  to  the  wall.  The  latter  aggravates  the  heat 
transfer  conditions  on  the  converging  wall  of  the  nozzle 
where  the  gaseous  heat  transfer  is  already  quite  severe. 

It  is  therefore  of  Interest  to  estimate  the  percentage 
of  impingement  of  the  particle  mass  flux  on  the  wall,  small  as 
it  may  be,  and  to  find  the  varioxis  parameters  of  significance. 

The  various  aspects  of  the  interaction  phenomena 
between  the  condensed  and  the  gaseous  phases  and  the  calculation 
of  the  thermodynamic  and  dynamic  properties  of  the  gas  and 
particles  are  discussed  in  detail  in  Reference  1.  The  results 
and  method  of  calculation  of  Reference  1  is  adopted  here  for 
investigating  the  inpingement  phenomena.  Accordingly,  the 


following  aBSuraptlonB  and  restrictions  made  in  Reference  1 
are  adopted: 

(a)  Small  pertur1:atlon  method.  The  convergeiKe 
of  the  solution  is  good  only  for  small  values 

of  the  small  disturbance  parameters  and  . 

(b)  Quasi  one'dlmenslonELl  assumption  of  gaseous  flow. 

(c)  Unifoim  distribution  of  particles  without 
interparticle  collisions. 

(d)  Stokes  drag  law  holds  for  the  micron  sl^^iC  particles. 

We  propose  to  use  several  different  approximte  metho«..s 
for  Investigating  the  inplngement  phenomena  to  reveal  the 
relative  significance  of  different  physical  factors.  Numerical 
calCTilatlons  of  the  differ«itial  equations  systems  governing 
the  impingement  will  also  be  performed  to  estimate  the 
overall  effect. 
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P..  THE  CONICAL  APPROXIMATION 


As  a  first  attempt  ve  shall  make  the  foUoving  assumptions, 
In  addition  to,  or  as  an  amplification  of  those  Bade  already  In 
Reference  1. 

(a)  The  gas  flow  Is  purely  conical,  l.e. ,  the  gas 
streamlines  are  rays  emanating  from  a  vertex 
downstream  of  the  throat. 

(b)  This  conical  flow  is  valid  from  the  entrance  of 

radius  r  to  the  point  N  where  the  tangent 
o 

to  the  geometric  throat  of  the  nozzle  meets  the 
extrapolated  conical  approximation  to  the 
nozzle  contour  (see  Figure  1). 

(c)  The  thermodynamic  and  flow  properties  of  both 
the  gas  and  the  particles  are  constant  adong 
circular  arcs  whose  center  is  in  the  vertex 
of  the  conical  surface. 

(d)  The  particles  will  possess  «  velocity  component 
perpendicular  to  the  conical  element,  l.e., 

a  tang«itlal  component.  There  will  be  no 
azimuthal  velocity  component  (at  least  on 
the  average).  The  radial  velocity  of  the  ^s 
and  the  particle  on  each  arc  defined  by  x 
(see  Figure  1)  are  those  determined  at  that  x 
by  the  one -dimensional  analysis  in  Reference  1. 

(e)  The  velfcity  of  the  particles  at  the  entrance  to 
the  nozzle  is  parallel  to  the  nozzle  axis,  i.e., 
the  initial  tangential  velocity  component  of 
the  particle  is 
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Seo  *  SrO-*«®o 


(see  Figure  1). 


We  shall  non-dimenslonallze  all  velocities  by  the 
speed  of  sound  in  the  combustion  chamber  ,  all  lengths  by 
the  length  of  the  converging  section  of  the  nozzle  L  ,  and 
time  by  a  reference  time 


thus  ve  get: 


r 


r  ■  — 
L 


z  » 


t 

K 


The  dlmensiorless  equations  of  motion  of  the  particles  in 
the  r  and  Q  directions  are  (see  Reference  2): 


(1) 


dz 


(2) 
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Centrifugal  tenns  In  eqns.  (l)  and  (2),  like  or 


were  neglected.  It  vas  checked  and  fbund  that  these  terms  did 
not  exceed  1^  of  the  terms  retained.  Also,  the  assumption  (c) 


vas  used  in  putting 


0  .  By  integrating  eqn.  (2)  and  defining 


L  de 


one  gets: 


1  “Z 

K  q  ^  •  -  e  dz 
0^90  r 


where 


o  L 


is  the  dimensionless  parameter  of  small  pr^rturbation  in  Reference  1, 

and  Is  the  Initial  tangential  vel  city  of  the  particle. 

We  shall  find  the  angle  9  from  which  the  pEurtlcIe  passing 

0 

through  the  point  N  originated.  The  area  percent  f  the 

annular  region  of  width  W  (Figure  1)  determined  by  9  ,  gives 

c  o 

a  close  upper  bound  of  the  percentage  of  pEirtlcle  flux  Ijqpinging 

on  the  converging  nozzle  wall: 

2nY  W  2W 

_  O  C  ,  C  /  1.  \ 


With  r  ■  r  -X  and 
o 
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.  -i  ^ 

"  K  dz 


together  with  our  asBumption  (e)  for  (3)  become b 


1  dx 

■  ^  °  Sr 


de 


S*tg0  • 

ro  o 


dx 


(5) 


By  Integrating  eqn.  (5)  fran  0  to  0  and  denoting  X  ■  r  -r. 

o  o  t 

and 

1  dx 

■  “  Sr 


X  ‘ 

/  - 


dx  »  I 


(6) 


one  getB  finally 


0-0  -  a  ‘tge  •  I  ■  0 


^ro 


(7) 


For  a  giwen  oonical  nozzle  and  known  valuea  of  K  and  q  , 

^  o  ^ro 

I  can  be  evaluated  and  0  1b  kiiown.  Eqn.  (T)  can  be  solved  for 

0^  either  graphically  or  by  trial  and  error. 

For  an  approxiiaate  soluulcn  of  0^  when  e  Ib  small 

take  tge  ■  0  thus 
o  o 


®0  “  ® 
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(8) 


also 


(9) 


W  ■  r  (8ln0-8in0  )  s  r  (0-e  ) 
c  o  o  o'  o 


using  eqns.  (8)  and  (9)  In  (^)  ve  get: 


P  s 
c 


2r^-© 


100 


(10) 


as  an  approxinate  expression  for  the  Impingement  percentage. 

Since  the  calculated  profiles  and  are 

taken  niunerlcally  as  functions  of  r  such  as  the  result  of 
the  computation  In  Reference  1^  it  is  possible  (a)  to  continue 
directly  the  nvuserlcal  program  and  e'valtiate  I  by  machine  confutation, 
or  (b)  q^  may  be  approximated  by  a  polyuomlal  of  some  degree  and 
then  Integrated  termvlse  analytically.  Eqn.  (10)  vLU  give  an 
estlnate  of  the  impingement  percentage.  This  approximate  method 
is  probably  adequate  vhen  a  substantial  portion  of  the  conrerging 
nozzle  is  approximately  conical. 
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3.  THE  N-S  METHOD 


In  this  section  ve  shall  adopt  the  streamlines  of  the 
gas  N  ■  const,  and  their  orthogonal  lines  S  ■  const,  as  the 
coordinate  system  (see  Figure  2).  This  N-S  coordinate  system 
offers  substantial  improvement  over  the  previous  conical  approxi¬ 
mation  near  the  entrance  and  near  the  throat.  We  shall  first 
keep  the  same  Initial  flow  conditions  entering  the  nozzle  for 
estimating  the  initial  tangential  velocity  of  the  particle  in 
order  to  focus  our  attention  on  the  role  of  the  shape  of  the  throat. 

The  one -dimensional  gas  and  particle  properties, 

either  thermodynamic  or  fluid  mechanical,  as  calculated  in 

Reference  1,  will  be  assumed  to  apply  along  S  •  const,  lines, 

l.e.  sdong  lines  perpendicular  to  the  gas  streamlines. 

Following  a  procedure  similar  to  that  adopted  in  the  conical 

approximation,  we  can  determine  the  path  and  coordinate  N 

(or  defined  in  eqn.  (15)  )  of  the  particle  at  the 

initial  section,  which  hits  the  vail  at  the  throat. 

N  or  t  determine  W_  (see  Figure  2)  which, under  the 
o  S 

assumption  of  uniform  particle  distribution  across  the  nozzle 
section,  gives  the  inpingement  percentage  as 


o  N.S. 


100  • 


Equations  (l)  and  ^2) 


S 

•  100  (11) 

Y 

o 

written  in  the  N-S  coordinates 


become: 
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(12) 


(13) 


Again  the  centrifugal  terms  in  eqn.  (12)  and  (13)  were  neglected 
for  the  same  reasons  given  in  eqn.  (l)  and  (2).  The  integral  of 
eqn.  (13)  with  the  use  of  the  assumption  of 


Sno  " 

gives 

1  dN 

Sn  "  K  dz  *  Sso  ^*^0  ^ 
o 


(14) 


Let  us  define  ^  as  the  stream  fiu«ction  of  the  gas  by: 


dt  -  P„(8)-q^(8)*c*d(N^) 


g 


(15) 


non-dimensionalization  with 

■  mass  flew  rate  of  gas 
g 

p*  •  gas  density  in  combustion  chaeber 
q*  ■  sonic  velocity  in  oonbustion  chanber 
L  »  le.ngt.h  of  converging  part  of  the  nozzle 


gives 


d\jf  ■  2ca 


Pg(8)q^(8)-KdH 


(16) 


.  9  - 


\rtiere 


p^q*L 


(0 


e 


with  i|r  ■  0  on  the  centerline,  and  i)r  ■  1  on  the  nozzle  vail. 

At  any  point  in  the  nozzle  defined  by  y'  and  S  ve  have 
the  relation: 


^  N(8,1)  ^  1.0 


(17) 


alao 


dt  ■  ^  *  T"  * 
^  dz 


(18) 


Using  eqna.  (17)^  (l6)  and  {ik)  in  (l8)  and  if  ve  vise  the  radiua 
y  instead  of  the  arc  length  N  with  less  tiian  1^  error  for 
cone  aisles  less  than  20*  (see  Figure  2)  ve  have  C  ■  k  in  (l6) 
one  finally  arrives  at: 

dt  -  2Ka-k^-q^tge^-p^(S)*q^(S)-N(S,l)*  if'i  ^dz  (19) 


Using  again  y(S,l)  instead  of  N(S,1)  ve  get  from  the  definition 
of  Tjf 

Pg(S)'q^(s)-*-y(s,l)^-a  •  1  (20) 

By  integrating  e<in.  (19)  from  to  1.0,  vith  the  help  of  eq^.  (20), 
one  gets: 
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(21) 


;^vii 

S8c‘«®0  ’  o 


Z,  Is  the  ti.De  needed  for  the  particle  vhich  starts  at  w  to 
1  o 

reach  the  throat.  Changing  the  variable  in  from  z  to  S 
through  the  definition  of  as  : 


O  0^8 


(22) 


ve  get  for  I, 


^throat  ^o  ® 

/  - -  dS 


(23) 


inust  usually  be  evaluated  numerically  because  of  the  appearance 
of  y(S.l)  ,  vhich  exp  xBases  the  knovn  nozzle  contour. 

Once  Is  evaluated,  eqn.  (21)  must  be  solved  for  . 

Pro.m  the  definition  of  ^  ve  have,  at  the  entrance  section: 


-  Q-p^(0)q^(0)*A^ 


where  A  ■  2icr^l  -CCS:?  )  is  the  spherical  entrance 


area 


(see  Figure  2)  which  results  in 


-  U 


(21^) 


CO80  ■  1-  , — 
o  h 


where 

h  -  apg(0)<i^(0)-2nr^  -  2  (~  )^. 

o 


(26) 


The  second  exprassion  for  h  employs  the  definition  of  a  . 
Eqn.  (21)  for  the  solution  of  becomes: 


so 


v/R 


^o  2 

(r ) 


1 


(2T) 


where 
has  to 


q  and  h 
^so 

be  evaluated 


are  known  for  a  certa.’?  problem  and 
from  eqn.  (23). 


I 


1 
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4.  TV/O-DIMENSIONAL  ENTRANCE  ANALYSIS 


Comparison  between  the  results  In  the  two  previous  sections 
serves  only  to  illustrate  the  function  of  tlrie  throat  shape .  Since 
the  same  approximtion  of  "conical"  entrance  condition  is  employed. 
With  a  major  part  of  the  impingement  taking  piace  in  the  entrance 
region  (to  be  shown  later).  We  have  to  investigate  carefully  the 
significance  of  the  geometry  of  the  entry  section.  To  facilitate 
comparison  witli  previous  results  we  shall  not,  however,  attempt 
a  strictly  axisymmetric  analysis.  We  shall  take  the  gas  velocity 
q  near  the  entrance  (see  Figure  3)  io  be  that  calculated  by  the 
use  of  the  one -dimensional  small  perturbation  analysis  in 
Reference  1.  Since  the  area  variation  near  the  entrance  is 
small,  this  is  a  very  good  assumption.  Adopting  the  cylindrical 
polar  coordinates  (x,y)  we  write  the  equations  of  social  (x)  and 
transverse  (y)  momentum  balance  of  the  particles  as: 


du 

P  ax 


r^Vp) 


(26) 


u 


dv 

dX 


—  (v  -V  ) 
K  ^  g  P^ 
o 


(29) 


In  order  to  solve  (28)  and  (2?)  for  the  particle  velocity  components 

u  and  V  ,  we  must  find  appropriate  expressions  for  the  gas 
P  P 

velocity  components  u  and  v  .  Since  the  angular  deflection 

g  g 

of  tne  gas  streamlines  is  small,  i.e.  p  is  small,  we  may  write 
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(30) 


u 

g 


q  ^oa^ 

g 


‘g 


is  known  from  tlie  one -dimensional  small  perturtation  analysis 
of  Reference  1.  Tne  transverse  velocity  component  of  gas 
(perpendicular  to  the  nozzle  eixis)  is: 


V 

g 


U  tgft  *  u 

g 


g  dx 


From  continuity  considerations  and  th.e  definition  of  the  stream 
function  ijf  we  get: 


y(x)  ^ 

Y(x) 

(see  Figure  3).  Using  eans .  (31)  in  (30),  we  get: 


(31) 


V 

g 


dY(x) 

dx 


(32) 


Y(x)  is  a  given  contour  function:  thus  U  and  v  are  known 

g  g 

functions  of  X  and  a  parameter  y  which  is  the  point  yhere  the 

o 

gaseous  element  originates.  With  the  help  of  eqns.  (3C>)  (32), 

the  Systran  of  two  ordinary  differential  equations  (26)  and  (29) 


serve  to  determine  U  and  v  as  functions  of  x  and  the 

P  P 

parameter  y^  .  To  solve  eqns.  (28)  and  (29)  by  power  series 
is  tedious  and  inaccurate. 


■* 
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They  are,  hence,  solved  simultaneously  on  a  computing 

machine  as  a  pe.rt  of  the  program  in  evaluating  q  from  the  siitail 

S 

pert.ur'bation  method  following  Reference  1.  Once  u  and  v  are 

P  P 

knovm,  they  can  be  integrated  in  the  same  program  to  give  the 
positior  and  direction  of  the  particle  at  any  station  x  as  a 
function  of  the  parameter  y^  as  : 


^  o 


o 


(33) 


where  C  is  a  constant  and  found  the  evalviation  of  the 
integral,  and 

V 

^  »  arc  tg  (~  )  (34) 

P 


The  percentage  of  Impingement  upstream  of  some  station 

X  is  found  by  determining  the  y^  of  the  pa.rticle  which  hits  the 

wail  at  that  X  .  This  y  will  be  denoted  by  y  .To  find 

o  po 

y  we  put  at  some  X  ,  y  ■  Y(x)  and  from  ecfi.  (33)  ve  have 
po  p 


(35) 


V  determines  the  width  of  impingement  region  and  again 

po  2D 

we  have  the  impingement  percentage  as 


•  100 


(36) 
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5 .  RESULTS  AND  COMPARISON  OF  THE  VARIOUS  fBTHODS 


From  the  studies  done  in  the  previous  sections  and  in 
Reference  1,  we  see  that  the  important  parameters  entering  the 
formulation  of  the  in?)inge?Tjent  are: 


2 


and  the  geometry  of  the  nozile.  Throu^out  the  ceLLculation 
two  inconsequential  assumptions  were  made  in  addition,  for  the 
sake  of  simplicity.  These  ar^: 

(a)  C  ■  C  gas  and  particle  specific  heats  are  equal, 

PP  Pg 

which  implies  that  ^  \  in  the  small  perturbation 

analysis . 

(b)  a  Kg  the  small  perturbation  pareuneter  for 
momentum  and  heat  transfer  are  equal.  All  small 
perturbations  were  calculated  to  the  first  order. 

In  order  to  follow  the  effect  of  each  individual 

parameter,  and  to  assess  the  validity  of  each  method,  many  cases 

were  calculated  with  results  summarized  in  Table  I.  Before 

analyzing  the  results,  let  us  first  get  some  idea  as  to  the 

lange  of  sizes  of  the  particles  for  \diich  our  analysis  may  be 

applicable.  Let  us  calculate  n  for  a  typical  value  of  k  ■  0,1 

3  “4 

L  ■  4",  fp  “  1*T2  grams/cm'^,  u  ■  6x10  poise,  U^  ■  10,000  cm/sec. 
Using  this  in 
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one  gets  t)  to  be  t)  ^  8  microns.  The  particle  sizes  corresponding 
to  other  values  of  in  a  nozzle  vith  the  same  L  are  given  .as 

Owe  range  of  Interest  is  hence  limited  to 

particles  of  only  several  microns  in  size  (less  than  13). 

A  comparison  between  cases  (i)  and  (2)  brings  out  the 
first  important  result.  The  results  obtained  in  case  (l)  using  the 
conical  method  and  those  in  case  (2)  using  the  N-3  method,  are  Ihe 
same.  Since  these  two  methoda  differ  only  in  the  throat  region,  we 
conclude  that  the  inpingement  is  insensitive  to  the  nozzle  shape 
near  the  throat.  This  fact  can  be  expiaiiaed  by  the  small 
relaxation  distance  of  the  exponential  decay  of  the  tangential 
velocity  expressed  by  eqns .  (3)  and  (l^+).  Because  of  this  rapid 
decay,  the  impingement  takes  place  mainly  in  the  entrance  section 
where  Z  is  small.  The  particle  motion  snon  follows  the  gas  motion 
upon  leaving  the  entrance.  Sizeable  changes  in  the  nozzle  shape 
near  the  throat  do  not  cheuige  appreciably  the  percentage  of 
impingement .  This  was  further  demonstrated  in  the  detailed 
results  of  the  numerical  conputations  of  cases  (3)  to  (20). 

Since  the  impingement  results  for  particles  of  sizes 
ranging  from  5-13  microns  depend  heavily  on  the  velocity  variations 
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of  the  gas  and  the  pairtlcles  tn  the  entrance  region,  It  is 
necessaxy  to  deteznlne  accurately  the  gas  and  the  particle 
motion  in  the  entrance  region.  The  method  of  area  perturbation 
as  is  adopted  in  Reference  1  is  inadequate  for  the  present  purpose 
Accordingly,  cases  (3)  to  (20)  were  calculated  for  a  prescribed 
nozzle  shape  (with  no  area  perturbation),  step  by  step  with 
the  help  of  an  IBM  7090  hi^  speed  computing  machine. 

Analysis  of  the  results  of  these  cases  shows: 

(a)  Moderate  changes  in  7  have  no  significant  effect 
on  the  impingement.  This  is  demonstrated  by  cases  (3) 
and  (4)  where  y  is  the  only  changing  parameter. 

The  results  of  these  two  cases  are  essentially  the  same. 

(b)  An  increase  of  the  particle  mass  flux  ratio  \ 
in  cases  (5)  to  (8)  slightly  decreases  the 
percentage  of  Impingement.  This  is  due  to  more 
heat  and  momentum  transfer  between  the  gaseous 

and  particle  phases,  which  affect  the  gaseous  velocity. 

(c)  When  dealing  with  a  prescribed  nozzle  geometry, 

the  most  important  parameter  affecting  the  amount  of 

impingement  is  the  particle  diameter  tj  The  effect 

of  TJ  can  be  detected  in  cases  (9)  to  (l4)  where  the 

only  varying  parsuneter  is  ■  Kg  •  As  was  shown 

before  for  a  fixed  nozzle  geometry  n  is  proportional 

to  .  The  results  of  cases  (9)  to  (1*^)  show  linear 

dependency  of  1  ,  which  determines  the  impingement, 

on  K  This  is  an  obvious  consequence  of  the 
o 

linearized  treatment  of  Reference  1.  Therefore, 
for  a  given  nozzle  the  impingement  percentage  is 
directly  proportional  to  the  square  of  the 
diameter  of  particles. 


(d)  Cases  (10)  to  (12)  and  (15)  to  (20)  give  the  effect 
of  the  increase  in  convergence  emgle  of  the  nozzle  from 
20® to  1*0®.  As  expected,  we  see  that  in5>lngeinent 
increases  with  an  Increase  of  convergence  angle. 

Althoui^  the  approximate  conical  method  is  relatively 
simple,  the  two-dimensional  method  yields  more  consistent  results. 
We  see  that  in  Nozzle  1  (cases  3*1^)  the  conical  results  are  almost 
always  higher  than  the  "tvD>*dlmensional"  results  >^le  in  Nozzles  2 
and  3  (cases  15-20)  the>  are  lower.  This  again  can  be  explained 
by  the  rapid  decay  of  tangential  velocity  in  the  "two-dimensional" 
analysis.  It  is  clear  that  the  latter  vhere  impingement  takes 
place  in  0.4 -0.5  of  the  nozzle  length  will  be  more  sensitive  to 
increase  in  convergence  angle  than  the  conical  method  where 
impingement  is  negligible  after  0.1  of  the  nozzle  length. 

All  the  results  summarized  in  Table  I  show  in  common 

that  the  overall  percentage  of  particle  mass  flux  that  impinge 

on  the  nozzle  wall  is  very  small.  In  the  range  of  interest  up 

to  K  »  0.2  ,  the  impingement  width  W  is  of  order  of  magnitude 
o 

of  the  displacement  thickness  of  the  turbulent  boundary  layer 
on  the  chanher  wall  at  the  nozzle  entrance.  In  the  boundary 
layer,  the  velocities  of  the  gas  and  the  particles  are  less  than 
those  given  by  the  quasi  cme -dimensional  results.  Ibe  relaxation 
distance  for  the  transverse  velocity  coi^nent  of  particle  will  be 
accordingly  smaller  and  the  particles  will  follow  more  closely  the 
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gas  streamline 8  with  less  impingement.  Hence,  the  present  results 
of  small  total  peirticle  impingement  is  actvially  an  upper  bound. 

In  this  case,  we  may  consider  our  results  as  an  upper  bound  for  the 
iniplngeiDent ,  end  conclude  that  in  this  si^e  range  impingement 
is  negligible. 

When  interpreting  cases  (17)  to  (20)  where  the  impingement 

is  somewhat  greater,  we  should  beeur  in  mind  that  tte  values  of  K 

o 

are  rather  large  and  the  vcdidity  of  the  linear  perturbation  results 
is  questionable.  These  cases  were  included  mainly  to  show 
a  qualitative  trend. 
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Figure  1.  The  conical  coordimtes 
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Figure  2.  The  N-S  ooordiiiates 
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Figxjre  k.  Conqparison  between  ccaical  and  two-dimenBlonal  methods 
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Figure  5»  Impingement  vs.  sgia.ll  parameter  K  . 
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Figure  6.  Nozzle  geometries 


